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We propose and analyze a new scheme to produce ultracold neutral plasmas deep in the strongly 
coupled regime. The method exploits the interaction blockade between cold atoms excited to high- 
lying Rydberg states and therefore does not require substantial extensions of current ultracold 
plasma experiments. Extensive simulations reveal a universal behavior of the resulting Coulomb 
coupling parameter, providing a direct connection between the physics of strongly correlated Ryd- 
berg gases and ultracold plasmas. The approach is shown to reduce currently accessible temperatures 
by more than an order of magnitude, which opens up a new regime for ultracold plasma research 
and cold ion-beam applications with readily available experimental techniques. 



Strongly coupled Coulomb systems occupy an exotic 
regime of plasma physics [1-3 where correlated dynam- 
ics due to strong interactions dominates random thermal 
motion of the charges. While plasmas generally consti- 
tute the most abundant state of matter [4 j, strongly cou- 
pled systems are rather scarce and typically require ex- 
tremely high densities that occur in such exotic settings 
as the interior of stars and giant gas planets [5 or inertial- 
confinement fusion experiments [6 . On the other hand, 
ultracold neutral plasmas (UNPs), produced by direct 
photoionization of laser-cooled atoms [ZH9], represent a 
promising alternative for creating and studying strongly 
coupled plasmas, which provides unique opportunities to 
probe dynamical phenomena [8 j. Their low tempera- 
tures yield strong coupling conditions at such small den- 
sities that typical evolution timescales are slowed down 
tremendously. This has enabled recent experiments to 
investigate collective excitations [lOj [11] , relaxation pro- 
cesses [I2j [13] as well as dynamical instabilities [14] at 
ultralow plasma temperatures. Such conditions also hold 
great promise for applications in nanotechnology, where 
UNPs are utilized to create high-brightness ion beams 
[T5H20] for fabricating and characterizing nanoscale ob- 
jects. 

Yet, these intriguing possibilities are currently lim- 
ited by intrinsic ion heating [2TH28] that takes place 
right after plasma creation. UNPs are produced from 
ultracold atoms whose initial temperatures suggest ex- 
tremely strong correlations, which originally raised high 
hopes for creating neutral plasmas with unprecedented 
coupling strengths [7 . However, this scheme generates 
the plasma far from equilibrium in an entirely disordered 
state, and the subsequent correlation build-up is accom- 
panied by substantial heating that pushes the plasma- 
ions just to the edge of the strong-coupling regime. Theo- 
retical work has suggested ways to suppress this disorder- 
induced heating (DIH) by introducing initial ion corre- 
lations, e.g., by pre-ordering atoms in an optical lattice 
[29] and by exploiting atomic correlations in a degenerate 
Fermi gas [21] or correlations due to Penning ionization 
of excited atoms [30 . These ideas set challenging ex- 
perimental demands such as perfectly filled lattices with 
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FIG. 1: (Color online) Schematics of the (a) conventional 
direct-ionization scheme for UNP production and (b) the two- 
step production proposed in this work. Starting from a ran- 
dom gas of ground state atoms (c), the latter creates a corre- 
lated arrangement of Rydberg atoms (d), whose subsequent 
ionization results in a correlated plasma (e). As shown in 
panel (f), these initial correlations substantially suppress DIH 
after plasma creation (solid line) compared to the conven- 
tional production scheme (dashed line). The depicted temper- 
ature evolution has been calculated for an initial plasma den- 
sity of p = 2 x 10 8 cm -3 and resonant excitation to Sr(100 1 5'o) 
Rydberg states (see text for laser parameters). 



~ 10 5 atoms or very high atomic densities, which, thus 
far, have precluded their experimental implementation. 

Here we describe a new UNP-production scheme that 
strongly suppresses DIH with currently available and well 
established experimental techniques. The approach is 
based on the strong correlations between Rydberg atoms 
that emerge when a cold gas of atoms is laser-excited to 
a high- lying electronic state. The frequency of the exci- 
tation laser turns out to permit continuous variation of 
the plasma coupling strength. For the first time, this will 
allow plasma temperatures and densities to be tuned in- 
dependently, enabling unprecedented studies of the phase 
diagram (31] of UNPs. Our calculations reveal a univer- 
sal relation between the plasma's coupling strength and 
the Rydberg excitation dynamics, suggesting that UNPs 
could provide a powerful probe of the correlation proper- 
ties of ultracold Rydberg gases. We demonstrate a max- 
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imum suppression of DIH by more than an order of mag- 
nitude, which opens up a new regime for UNP research 
and significantly boosts the brightness and resolution of 
UNP-based ion beam technology [15]. 

Figure [I] illustrates the basic idea in comparison to 
the conventional approach for UNP creation. In the lat- 
ter case, direct photoionization of cold atoms [Figjlja)] 
produces a disordered configuration of ions, whose sub- 
sequent rearrangement establishes inter-particle correla- 
tions, and thereby decreases the potential energy of the 
plasma. Since the total energy is conserved, this leads to 
a rapid rise of the temperature to TJ ~ 1 K, as shown in 
FigjTJf). The corresponding degree of equilibrium cor- 
relations can be characterized by the so-called Coulomb 
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coupling parameter T = q ^ bT , where e is the electron 

charge, &b is the Boltzmann constant, and a = (|7rp)~i 
is the Wigner-Seitz radius for a plasma of density p. A 
plasma is termed strongly coupled, when the average po- 
tential energy ~ e 2 /a of the charges exceeds their thermal 
energy ~ k^T, i.e. when T > 1. Starting from a random 
distribution of almost stationary ions, the subsequent 
equilibration and DIH establishes a Coulomb coupling 
parameter of Y w 2 irrespective of the initial ion density 
and atom temperature [HJ [25] [28] [cf. Fig{TJf)]. This 
strong heating can be reduced with a two-step ioniza- 
tion scheme [see Figfljb)], where the ground state atoms 
are first excited to high-lying Rydberg states and sub- 
sequently ionized. In the first excitation step we take 
advantage of the so called Rydberg blockade [32], which 
prevents simultaneous excitation of nearby atoms due to 
the enormous van der Waals interaction between Ryd- 
berg atoms [33]. As illustrated in Figjljd), this gives 
rise to strong Rydberg-Rydberg atom correlations (see 
also Fig{2|. Their subsequent ionization thus produces a 
pre-correlated plasma [Figjlje)] that evolves to an equi- 
librium state with a much higher Coulomb coupling pa- 
rameter [Figjljf )] . Indeed the final plasma correlation 
function, shown in Figj2] closely resembles that of the 
prepared Rydberg gas but also shows pronounced oscilla- 
tions characteristic of a strongly coupled Coulomb liquid. 

The Rydberg blockade effect has been the basis of 
recent experimental breakthroughs for applications in 
quantum information [34 36] and nonlinear quantum op- 
tics [37-42 , and opened the door for studies of quan- 
tum many-body phenomena [43H48] . While most of 
these experiments utilize alkaline atoms, we focus here 
on cold Strontium atoms whose secondary valence elec- 
tron provides a viable way to probe the temperature of 
the created plasma ions via optical absorption and fluo- 
rescence measurements [23j 26, 49 . Following recent ex- 
periments on Sr-Rydberg gases [50-52 , we consider res- 
onant two-photon excitation from the Sr(5 1 5'o) ground 
state to a Sr(n 1 S'o) Rydberg state via the intermediate 
Sr(5 1 Pi) state by two excitation lasers with Rabi fre- 
quencies Vt\ and 0,2 and a two-photon detuning A [see 




FIG. 2: (Color online) (a) Schematic level diagram for two- 
photon Rydberg excitation, (b) Correlation function g(r) for 
the Rydberg ensemble and the resulting plasma after relax- 
ation for the same parameters as in FigJTJf ) . The vertical line 
marks the blockade radius according to eq.Q. (c) Correla- 
tion function of the Rydberg ensemble for two different laser 
detunings A. 



Fig{2|a)]. The single-atom steady state excitation frac- 
tion P ry d = + needs to be high in order 
to maximize interaction effects. We, therefore, choose 
fti = 25MHz and ^ 2 = 0.5MHz, which gives P ryd « 1 
and can be realized in current experiments [50j [51] . Due 
to the large decay rate 7 = 32MHz of the intermediate 
state, the resulting many-body excitation dynamics can 
be described in terms of transition rates [53] [54]. The cor- 
responding set of rate equations can be solved efficiently 
via classical Monte Carlo simulations that account for 
the Rydberg-Rydberg atom interaction via an effective 
detuning [53] [5l] 

- C 6 



A 7 = A ■ 



(i) 



of the ith atom at position r^, where the sum runs over all 
atoms that are in the Rydberg state. This approach has 
been applied to study cold Rydberg gases [53 -55 and lat- 
tices [53 , 55 -57] and was shown to yield good agreement 
with quantum calculations [54, 56 . For large principal 
quantum numbers, n, the van der Waals coefficient Cq of 
Sr(n 1 6'o) Rydberg states [58, 59 takes on enormous val- 
ues such that a single Rydberg atom can inhibit the exci- 
tation of up to a few 10 3 atoms. This implies a large num- 
ber of ground state atoms required to obtain sufficiently 
many Rydberg excitations. Such sizable ensembles can 
be described by the Monte Carlo procedure, which al- 
lows us to treat up to several 10 7 atoms. Typical simu- 
lation results are illustrated in Figj2] where we show the 
Rydberg-Rydberg atom correlation function after laser 
excitation. At short distances, simultaneous excitation is 
completely suppressed by the interaction blockade [32] . 
The corresponding blockade radius 
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can be estimated by equating the respective interaction 
energy with the linewidth Az/ [85] of the Rydberg ex- 
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FIG. 3: (Color online) (a) Considered shape of the field pulse 
for Rydberg atom ionization with r = 2tt/5uj. (b) Ioniza- 
tion fraction, /i on , and (c) electron excess energy, E e , for 
(27TCJ)- 1 = 100 ps, n = 70 and different values of Fq. 



namics shown in FigjTJf), which demonstrates a more 
than tenfold enhancement of the Coulomb coupling pa- 
rameter. This places the produced plasma deep into the 
strongly coupled regime with strong liquid like correla- 
tions [see Fig{2|a)]. 

The described production scheme involves several pa- 
rameters that can be widely tuned experimentally. To 
elucidate their effects we have determined the achiev- 
able coupling strengths T for various combinations of the 
ground state atom density po? the addressed Rydberg 
level n and the detuning A of the Rydberg excitation 
laser. For fixed t exc = 0.5/is and resonant excitation, the 
Coulomb coupling parameter T increases with both the 
atomic density po and the principal quantum number n 
of Rydberg state [FigQa)]. This behavior can be quali- 
tatively understood from the increasing degree of corre- 
lations in cold Rydberg gases with po and the strength 
of the van der Waals interactions. We can quantify this 
relation by defining the number of blockaded atoms 



citation. This simple relation yields a good descrip- 
tion of the numerical results in Figj2]and generally gives 
r*b ~ 1.2 a, which is close to the typical correlation length 
of a strongly coupled plasma. 

The conversion of these atomic correlations into ionic 
ones, requires efficient ionization of the Rydberg states 
while leaving the ground state atoms undisturbed. To 
this end, we consider ionization by a short electric field 
pulse, which we describe by classical trajectory Monte- 
Carlo (CTMC) simulations [60]. We assume a simple 
pulse shape F = F cos(u;t + <p) exp(— t 2 /r 2 ), where Fq 
denotes the field amplitude and cj, <\> and r parametrize 
the form of the pulse [cf. Fig{3^a)]. The ionization pro- 
cesses has to fulfill two major requirements: (i) ionization 
needs to be efficient to maximize the transfer of the Ry- 
dberg atom ordering into ionic correlations and (ii) the 
excess energy E e of the ionized electrons should be suffi- 
ciently low so that they remain trapped in the collective 
ionic space charge potential to form a neutral plasma 
state [7 j. As shown in Figj3^b) both conditions are opti- 
mized around <fi ~ 7r/2, where we find a low electron en- 
ergy - 100K and nearly 100% ionization of Sr(70%) Ry- 
dberg states. This energy is low enough to enable plasma 
formation [7 but sufficiently high to prevent rapid re- 
combination after plasma creation [61-65 . Such ideal 
conditions can also be obtained with other pulse shapes, 
such as, e.g., half-cycle pulses as shown theoretically and 
experimentally in [66 . 

Finally, we performed molecular dynamics simula- 
tions [67] to investigate the subsequent relaxation of 
the plasma, which is described as a one-component sys- 
tem [86 of initially stationary ions, whose initial po- 
sitions correspond to those of the Rydberg atoms ob- 
tained from the Monte-Carlo simulations of the excita- 
tion step. Such simulations yield the temperature dy- 
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which corresponds to the total number of atoms within 
one blockade sphere of radius ?v Indeed all data points 
collapse on a single universal curve as a function of iVb 
[see Figj4^b)]. In the strong blockade limit V shows a 
weak logarithmic growth with N\> such that further in- 
crease of the ground state atom density does not yield a 
substantial enhancement. 

However, the detuning of the Rydberg excitation laser 
provides another control parameter to tune the Coulomb 
coupling strength. This is demonstrated in Figj5] where 
we show T as a function of A for different ground state 
densities and Rydberg states. Apparently the univer- 
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FIG. 4: (Color online) (a) Final Coulomb coupling parameter 
T as a function of the ground state atom density po for t exc = 
0.5 ps and different Rydberg levels n. (b) Same data, but as 
a function of the number iVb [eq.Q] of blocked atoms. The 
inset demonstrates the logarithmic scaling (dashed line) of T. 



4 





n 


p [cm 3 ] 




• 100 


10 12 


5160 


▲ 100 


2.9xlO u 


1500 


V 90 


5.3x1011 


1500 


■ 70 


7.5x1011 


500 



20 -10 
A/2tt [MHz] 



FIG. 5: (Color online) Coulomb coupling parameter T as a 
function of laser detuning A for t e xc = 0.5 /is and different 
combinations of n and po. The horizontal dotted line shows 
the Coulomb coupling parameter for direct ionization. 



sal scaling also holds off resonance, such that the final 
result is determined by A and only. For the at- 
tractive interactions between Sr(n 1 5'o) atoms [58] [59] a 
finite laser detuning gives rise to interaction-induced res- 
onances at the red side of the excitation line (A < 0), 
where Rydberg excitation is enhanced for configurations 
where A^ « 0. This additional distance-selectivity can 
produce a more pronounced peak structure of the corre- 
lation function [68] and tends to steepen its drop around 
rb [see Figj^b)]. As a result the Coulomb coupling pa- 
rameter can be increased by tuning to the red side of the 
atomic resonance. For large TVb it assumes a maximum 
around A = —7/2 and drops at either side of the atomic 
line. Well off resonance, DIH is even enhanced, such that 
T can be tuned from the weak-to-moderate coupling re- 
gion deep into the strong coupling regime with maximal 
Coulomb coupling parameters around T ~ 35. 

The proposed scheme, therefore, represents a viable 
way to not only tremendously boost currently achiev- 
able coupling strength but also to independently tune the 
plasma parameters of UNPs, which would otherwise be 
constrained to V = 2 oc p 1 / 3 T. The demonstrated tenfold 
cooling directly translates into an enhanced brightness of 
UNP-based ion beams, which enhances the spatial beam 
resolution [15] for nanotechnology applications. 

We have focussed on electric-field ionization of the Ry- 
dberg gas, but other approaches for plasma-conversion 
should also be possible, such as long- wavelength pho- 
toionization or electron impact ionization. As pointed 
out recently [69] , the latter occurs naturally in dense Ry- 
dberg gases, driving a spontaneous avalanche-like evo- 
lution to a plasma [7QH74] that could inherit residual 
Rydberg atom correlations. Spontaneous plasma forma- 
tion, however, takes place on long microsecond timescales 
[69j [70] and involves significant atomic motion [75] and 
loss of highly-excited atoms to low-lying states [TU [72] 
— all of which deteriorate the cooling effect described 
here. However, electron impact ionization can be driven 
with less deleterious effect by combining Rydberg exci- 
tation with the conventional method for plasma creation 



[7] in a double-pulse scheme. Here, the first laser-pulse is 
tuned below threshold to excite strongly correlated Ry- 
dberg atoms while leaving the majority of atoms in the 
ground state. The second pulse, tuned above thresh- 
old, ionizes a fraction of the ground state atoms to pro- 
duce a low-density plasma with energetic electrons that 
subsequently convert the Rydberg gas into a correlated 
plasma. Such double-pulse sequences have already been 
implemented for pump-probe experiments in Rydberg 
gases [44] [76] , to excite Rydberg atoms in a UNP [77] [78] , 
and to convert a Rydberg gas to a plasma [51] . The final 
degree of plasma correlations will depend on the parame- 
ters of the seed-plasma, which raises interesting questions 
for future studies. 

Several extensions of the proposed scheme appear 
worth pursuing. First, a correlated plasma with T < 
lOOmK seems well suited for subsequent laser-cooling on 
the Sr + (5Si/2) Sr + (5Pi/2) transition, which thus far 
was hampered by DIH out of the Doppler range. Sec- 
ond, long-range dipolar Rydberg interactions, induced by 
resonant pair-state couplings [79-81 or external electric 
fields [82 , may give rise to atomic correlations that are 
closer to those of long-range interacting Coulomb sys- 
tems and, therefore, may further enhance the coupling 
strengths of the resulting plasma. Here, higher order in- 
teractions give rise to an interesting, correlated excitation 
dynamics [8TJ [83l EI] that could be probed by monitor- 
ing T of the corresponding UNP. Finally, Rydberg excita- 
tion by standing-wave light fields is expected to increase 
atomic correlations substantially. If the lattice constant 
matches the blockade radius, we expect strong excitation 
ordering that could push T to values not far below the 
crystalline-plasma regime (T & 174 pQ). 
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